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Abstract 

The bulk of the hot and dense matter created at RHIC behaves like an almost 
ideal fluid. I present the evidence for this and also discuss what we can learn 
about the transport properties of the quark-gluon plasma (QGP) from the grad¬ 
ual breakdown of ideal fluid dynamic behavior at large transverse momenta, 
lower beam energies, larger impact parameters, and forward rapidities. 

1 The QCD Equation of State and ideal fluid dynamics 

With relativistic heavy-ion collisions one explores the phase diagram of strongly interacting bulk matter 
in the regime of high energy density and temperature. Lattice QCD (LQCD) tells us (T) that for zero net 
baryon density QCD matter undergoes a phase transition at T cr = 173 ± 15MeV from a color-confined 
hadron resonance gas (HG) to a color-deconfined quark-gluon plasma (QGP). The critical energy density 
e cr ~ 0.7 GeV/fm 3 £Q| corresponds roughly to that in the center of a proton. At the phase transition, the 
normalized energy density e/T 4 rises rapidly by about an order of magnitude over a narrow temperature 
interval AT < 15 — 20MeV, whereas the pressure p/T 4 (which is proportional to the grand canonical 
thermodynamic potential) is continuous and rises more gradually (Fig.|T|). Both seem to saturate at about 
80-85% of the Stefan-Boltzmann value for an ideal gas of noninteracting quarks and gluons, the energy 
density more quickly (at about 1.2 T cr ), the pressure more slowly. Above about 2 T cr , the lattice data 
follow the Equation of State of an ideal gas of massless particles, e = 3 p. 

For many years this observation has been interpreted as lattice QCD support for the hypothesis of 
a weakly interacting, perturbative QGP. The recent RHIC data taught us that this interpretation was quite 
wrong. The first part of the title of this talk, which was lifted from a coffee mug nowadays distributed by 
Brookhaven National Laboratory to their guests, alludes to this exciting discovery. 

It was recognized over 3 decades ago (see review El) that information about the EOS of strongly 




Fig. 1: The normalized energy density e/T 4 (left) and pressure p/T 4 (right) from lattice QCD JT| for 0, 2 and 3 
light quark flavors, as well as for 2 light + 1 heavier (strange) quark flavors. Horizontal arrows on the right indicate 
the corresponding Stefan-Boltzmann values for a non-interacting quark-gluon gas. 

interacting matter can be extracted by studying the collective dynamics of relativistic heavy-ion colli¬ 
sions. This connection is particularly direct in the framework of ideal fluid dynamics which becomes 
applicable if the matter formed in the collision approaches local thermal equilibrium. The latter requires 
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sufficiently strong interactions in the medium that local relaxation time scales are shorter than the macro¬ 
scopic evolution time scale. In this limit the local conservation laws for the baryon number, energy and 
momentum currents, d fl j^(x) = 0 and d l ,T ,w = 0, can be rewritten as the relativistc Euler equations for 
ideal fluid motion: 


riB = -riB(d-u), e = - 

ie + p) (d-u), 

(1) 

*<* - V “ V - C * , V'Unf 

-)■ 

(2) 

e + p 1 + cj \ 
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The dot denotes the time derivative in the local fluid rest frame (/ = u ■ df) and V /( the gradient in the 
directions tranverse to the fluid 4-velocity vl‘. The first line describes the dilution of baryon and energy 
density due to the local expansion rate d ■ u, which itself is driven according to Q by the pressure or 
energy density gradients providing the fluid acceleration. The absolute value of the energy density e is 
locally irrelevant: the initial maximal energy density eo only matters by setting the overall time scale 
between the beginning of hydrodynamic expansion and final decoupling, thereby controlling how much 
flow can develop globally. The details of the flow pattern are thus entirely controlled by the temperature 
dependent speed of sound c 2 = j^. 




Fig. 2: The square of the speed of sound cj. from lattice QCD above T cr ffl (left) and from models above and below 
T cr 0 (right). 

According to the LQCD data, the latter is c 2 « 2 for T > 2 T cr , then drops steeply near T ~ T cr 
to values near c 2 « (the “softest point”, see Fig. |2j left panel), before rising again in the hadron 
resonance gas phase to c 2 « 0.15 0 (Fig.0 right panel). A key goal of flow studies in relativistic heavy 
ion collisions is to find traces of this “softest point” in the data. 

2 “Flavors” of transverse flow in heavy ion collisions 

Experimentally one studies flow by analyzing the transverse momentum spectra of the emitted hadrons. 
In central (b = 0) collisions between spherical nuclei, the flow is azimuthally symmetric about the beam 
axis. This “radial flow” integrates over the entire pressure history of the collision from initial thermal- 
ization to final decoupling (“freeze-out”), due to persistent pressure gradients. In noncentral (b / 0) col¬ 
lisions, or central collisions between deformed nuclei such as uranium 0, the nuclear reaction zone 
is spatially deformed, and anisotropies of the transverse pressure gradients result in transverse flow 
anisotropies. These can be quantified by Fourier expanding the measured final momentum spectrum 
dN/(dy p±dp±d(j)p) with respect to the azimuthal angle fl p . For collisions between equal nuclei, the 
first non-vanishing Fourier coefficient at midrapidity is the elliptic flow t’a ip ±, b). Since V 2 is driven by 
pressure anisotropies and the spatial deformation of the reaction zone creating such anisotropies quickly 
decreases as time proceeds, the elliptic flow is sensitive to the EOS only during the early expansion stage 
llhl (the first ~ 5 fm/c in semicentral Au+Au collisions 0 ), until the spatial deformation has disappeared. 

Depending on the initial energy density (i.e. beam energy), the hot expanding fireball spends 
this crucial time either entirely in the QGP phase, or mostly near the quark-hadron phase transition, or 
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predominantly in the hadron resonance gas phase (7J, thereby probing different effective values of the 
sound speed c s . To the extent that ideal fluid dynamics is valid in all these cases, an excitation function 
of the elliptic flow V 2 should thus allow to map the temperature dependence of the speed of sound and 
identify the quark-hadron phase transition, via a minimum in the function (\fs) Q This will be 
further discussed below (see Section l3~2l and Fig. |TJi. 

3 Model parameters and predictive power of hydrodynamics 

The hydrodynamic model requires initial conditions at the earliest time at which the assumption of local 
thermal equilibrium is applicable, and a “freeze-out prescription" at the end when the system becomes 
too dilute to maintain local thermal equilibrium. Both are described in detail elsewhere 0. Different 
approaches to freeze-out invoke either the Cooper-Frye algorithm |9) (used by us), in which chemical 
freeze-out of the hadron abundances at T cr fTO) must be implemented by hand by introducing non¬ 
equilibrium chemical potentials below T cr I I1111211X31 . or a hybrid approach [.14j,T5] that switches from 
a hydrodynamic description to a microscopic hadron cascade at the quark-hadron transition, letting the 
cascade handle the chemical and thermal freeze-out kinetics. While the radial flow patterns from the two 
freeze-out procedures don’t differ much, the elliptic flow can be quite different if the spatial deformation 
of the fireball is still significant during the hadronic stage of the expansion, as I will discuss in Sec. 15.21 

We have solved the relativistic equations for ideal hydrodynamics under the simplifying assump¬ 
tion of boost-invariant longitudinal expansion (see dm for details). This is adequate near midra¬ 
pidity (the region which most RF1IC experiments cover best), but not sufficient to describe the rapidity 
distribution of emitted hadrons and of their transverse flow pattern which require a (3+l)-dimensional 
hydrodynamic code such as the one by Hirano HU . 

The initial and final conditions for the hydrodynamic evolution are fixed by fitting the pion and 
proton spectra at midrapidity in central collisions; additionally, we use the centrality dependence of the 
total charged multiplicity dN^/dy. I stress that this is the only information used from bf 0 collisions, 
and it is necessary to fix the ratio of soft to hard collision processes in the initial entropy production. 
An upper limit for the initial thermalization time tq < 0.6 fm/c, the initial entropy density so = 110 fm 3 
in the fireball center (corresponding to an initial peak energy density eo ~ 30 GeV/fm 3 and central fire¬ 
ball temperature To ~ 360 MeV ~ 2 T cr ), the baryon to entropy ratio, and the freeze-out energy den¬ 
sity edec = 0.075 GeV/fm 3 are all fixed from the 6 = 0 pion and proton spectra (these numbers refer to 
200 A GeV Au+Au collisions at RHIC fi~3l ). The initial entropy density profile is calculated for all b 
from the collision geometry, using the Glauber model with soft/hard ratio as fixed above. Except for pi- 
ons and protons, all other hadron spectra in b = 0 collisions and all spectra for b / 0 collisions (including 
all flow anisotropies such as v 2 which vanish at b = 0) are then parameter-free predictions of the model. 
Note that all calculated hadron spectra include feeddown from decays of unstable hadron resonances. 

4 Successes of ideal fluid dynamics at RHIC 
4.1 Hadron momentum spectra and radial flow 

Figure |3] shows the single particle pj_-spectra for pions, kaons and antiprotons (left panel) as well as G 
baryons (right panel) measured in Au+Au collisions at RHIC, together with hydrodynamical results 1131 . 
In order to illustrate the effect of additional radial flow generated in the late hadronic stage below T cr , 
two sets of curves are shown: the lower (blue) bands correspond to kinetic decoupling at T cr = 165 MeV, 
whereas the upper (red) bands assume decoupling at T ( \ ec = 100 MeV. The width of the bands indicates 
the sensitivity of the calculated spectra to an initial transverse flow of the fireball already at the time of 
thermalization fT3l In the hydrodynamic simulation it takes about 9-10 fm/c until most of the fireball 
becomes sufficiently dilute to convert to hadronic matter and another 7-8 fm/c to fully decouple Q. Fig¬ 
ure 0 shows that by the time of hadronization the dynamics has not yet generated enough radial flow to 
reproduce the measured p and G spectra; these heavy hadrons, which are particularly sensitive to radial 
flow, require the additional collective “push” created by resonant quasi-elastic interactions during the 
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at RHIC m- The curves show hydrodynamical calculations O (see text). 

fairly long-lived hadronic rescattering stage. The flattening of the p spectra by radial flow provides a nat¬ 
ural explanation for the (initially puzzling) experimental observation that for p±>2 GeV/c antiprotons 
become more abundant than pions ED- 

As shown elsewhere (see Fig. 1 in 1181 ). the model describes these and all other hadron spectra 
not only in central, but also in peripheral collisions, up to impact parameters of about 10 fm, and with 
similar quality. No additional parameters enter at non-zero impact parameter. 

4.2 Elliptic flow 

Figure |4] shows the predictions for the elliptic flow coefficient v-i from Au+Au collisions at RHIC, to¬ 
gether with the data lTT9l [20l . For impact parameters h< 7fm (n c h/n max >0.5) and transverse momenta 




Transverse Momentum p T (GeV/c) 


Fig. 4: Left: p -averaged elliptic flow for all charged hadrons from 1 30 A GeV Au+Au collisions, as a function 
of collision centrality (n c h is the charged multiplicity at y = 0). The curves are hydrodynamic calculations with 
different choices for the initial energy density profile (see 12T1 ). Right: Differential elliptic flow V 2 (p±) for identi¬ 
fied hadrons from minimum bias Au+Au collisions at 200 A GeV 1 19il20lf22] |. together with hydrodynamic curves 
from l23l . 

p l < 1.5 GeV/c the data are seen (left panel of Fig. EJ to exhaust the upper limit for +2 obtained from 
the hydrodynamic calculations. For larger impact parameters b > 7fm the p L -averaged elliptic flow +2 
increasingly lags behind the hydrodynamic prediction; this will be discussed in detail in Sec. 15.21 As a 
function of p_ (right panel of Fig. Eli the elliptic flow of all hadrons measured so far is very well described 
by hydrodynamics, for p± < 1.5 — 2 GeV/c. In particular the hydrodynamically predicted mass splitting 
of V 2 at low p± is perfectly reproduced by the data. This mass splitting depends on the EOS G2), and 
the EOS including a quark-hadron phase transition used here describes the data better than one without 
phase transition (see Fig. 2 in Ideal fluid dynamics with a QGP EOS thus gives an excellent and 
very detailed description of all hadron spectra below p± = 1.5 GeV/c. Since this p± -range includes more 
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than 99% of all produced hadrons, it is fair to say that the bulk of the fireball matter formed in Ait+Au 
collisions at RHIC behaves very much like a perfect fluid. 

4.3 Final source eccentricity in coordinate space 

While spectra and elliptic flow reflect the momentum structure of the hadron emitting source, its spatial 
deformation can be tested with 2-pion Hanbury Brown-Twiss (HBT) correlations measured as a function 
of the azimuthal emission angle [24 25 j|. Even though the initial spatial deformation of the reaction zone 
in non-central Au+Au collisions at RHIC is finally completely gone, many pions are already emitted 
from earlier times when the spatial deformation is still significant. For Au+Au collisions at b = 7 fm, 
the spatial eccentricity of the time-integrated hydrodynamic pion emission function is (e x ) = 0.14 (or 
56% of its initial value e x {jf) = 0.25) l24l . Using azimuthally sensitive pion HBT measurements, the 
STAR Collaboration has measured in the corresponding impact parameter bin ll26l (e x ) = 0.11, ±0.035 
(or 45 ± 15% of the initial deformation). This can be counted as another success for hydrodynamics. 

5 Viscous effects at RHIC 
5.1 QGP viscosity 

As evident in the right panel of Fig.|4] the hydrodynamic prediction for '(’2 (p±) gradually breaks down 
above p± > 1.5 GeV/c for mesons and above p± > 2.2 GeV/c for baryons. The empirical fact l22l that 
both the values, where this break from hydrodynamics sets in, and the saturation values for t ?2 at high 
p± for baryons and mesons are always (;. e. for all collision centralities!) related by 3 :2 (i.e. by their 
ratios of valence quark numbers), independent of their masses, tells its own interesting story (see e.g. 
Erl '): It strongly suggests that in this p± region hadrons are formed by coalescence of color-deconfined 
quarks, and that the elliptic flow is of partonic origin (i.e. generated before hadronization), with a p_- 
shape that follows hydrodynamics at low pj_ up to about 750 MeV and then gradually breaks away E3- 

Since V 2 {p±) is a measure for the relatively small differences between the in-plane and out-of- 
plane slopes of the p± spectra, it is more sensitive to deviations from ideal fluid dynamic behaviour than 
the angle-averaged slopes. Two model studies [28,29 ] showed that 02 reacts particularly strongly to shear 
viscosity. As the mean free path of the plasma constituents (and thus the fluid’s viscosity) goes to zero, 
V 2 approaches the ideal fluid limit from below f30l (see Fig.[5j0- At higher transverse momenta it does 
so more slowly than at low p± approaching a constant saturation value at high p±. The increasing 
deviation from the ideal fluid limit for growing p± is qualitatively consistent with the expected influence 




pj_ [GeV] 

Fig. 5: Left: Elliptic flow from a parton cascade |!30j|, compared with STAR data, for different parton-parton 
scattering cross sections. Larger cross sections lead to smaller mean free paths. Right: Perturbative effects of shear 
viscosity on the elliptic flow V 2 (p±) 1291 (see text). 

of shear viscosity: Teaney f29l showed that the lowest order viscous corrections to the local equilibrium 
distribution increase quadratically with p± so that V 2 remains increasingly below the ideal fluid limit as 
72 j grows (see Fig. Eb). From Fig. Eb Teaney concluded that at RHIC the normalized sound attenuation 
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length = rpj-z 2 (where r/ is the shear viscosity, T the temperature and s the entropy density) cannot 
be much larger than about 0.1. This puts a stringent limit on the dimensionless ratio r//s, bringing it 
close to the recently conjectured absolute lower limit for the viscosity of p/s = ?i/(47r) |3TJ- This would 
make the quark-gluon plasma the most ideal fluid ever observed HQ. 

These arguments show that deviations from ideal fluid dynamics at high p± must be expected, and 
that they can be large even for fluids with very low viscosity. At which p± non-ideal effects begin to 
become visible in V 2 (p±) can be taken as a measure for the fluid’s viscosity. To answer the quantitative 
question what the RHIC data on partonic elliptic flow and its increasing deviation from ideal fluid be¬ 
haviour above p± > 750 MeV/c imply for the value of the QGP shear viscosity p requires a numerical 
algorithm for solving viscous relativistic hydrodynamics, see f32l . 

5.2 Viscosity of the hadron resonance gas 

Ideal fluid dynamics also fails to describe the elliptic flow V 2 in more peripheral Au+Au collisions at 
RHIC and in central and peripheral collisions at lower energies (see Fig. St), as well as at forward 
rapidities in minimum bias Au+Au collisions at RHIC EH- Whereas hydrodynamics predicts a non¬ 
monotonic beam energy dependence of V 2 (Fig.[7Ji Q), with largest values at upper AGS and lower SPS 
energies, somewhat lower values at RHIC and again larger values at the LHC, the data seem to increase 
monotonically with y/s. 


OJ 




Fig. 6: Left: Scaled elliptic flow (e denotes the initial spatial eccentricity) vs. the charged multiplicity per 
unit initial transverse overlap area S E3- Right: Elliptic flow i >2 for minimum bias Au+Au collisions at various 
collision energies, parametrized by the final charged multiplicity density at midrapidity ca 

Such a monotonic rise is consistent with “hybrid” calculations by Teaney (31 (Fig. Hr) where the 
fireball undergoes ideal fluid dynamic evolution only while in the QGP, followed by hadronic kinetic 
evolution using RQMD after hadronization. Figure [dji shows several curves corresponding to different 
equations of state during the hydrodynamic evolution (see I I51 ). with LH8 being closest to the lattice 
data. The difference between the points labelled LH8 and the hydrodynamic values at the top of the figure 
is due to the different evolution during the late hadronic stage. Obviously, at lower collision energies and 
for impact parameters b ~ 7 fm (simulating minimum bias collisions), ideal fluid dynamics continues 
to build additional elliptic flow during the hadronic stage, but RQMD does not. The initial energy 
densities are smaller than at RHIC and the fireball does not spend enough time in the QGP phase for the 
spatial eccentricity s to fully relax before entering the hadron resonance gas phase. Anisotropic pressure 
gradients thus still exist in the hadronic phase, and ideal fluid dynamics reacts to them according to the 
stiffness of the hadron resonance gas EOS (p ~ 0.15e). Teaney’s calculations fl5l show that RQMD 
responds to these remaining anisotropies much more weakly, building very little if any additional elliptic 
flow during the hadronic phase. The hadron resonance gas is a highly viscous medium, unable to maintain 
local thermal equilibrium. The failure of the hydrodynamic model in situations where the initial energy 
density is less than about 10GeV/fm 3 l35l is therefore likely not the result of viscous effects in the early 
QGP fluid, but rather caused by the highly viscous late hadronic stage which is unable to efficiently 
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respond to any remaining spatial deformation. This is supported by a compilation of midrapity data on 
single particle spectra and elliptic flow as a function for p± for pions and protons by the PHENIX Colla¬ 
boration (Fig. 20 in El) which shows that (i) no purely ideal fluid dynamic model can describe all the 
data, (ii) differences among theories and between theoretical predictions and the experimental data can 
all be traced to different ways of describing the hadronic stage of the expansion, and (iii) the only model 
which describes all data simultaneously is the hybrid hydro+cascade approach by Teaney et al. El- 

Similar arguments hold at forward rapidities at RHIC <351 where the initial energy densities are 
also significantly smaller than at midrapidity while the initial spatial eccentricities are similar. We re¬ 
cently performed an analysis of the rapidity dependence of the charged hadron elliptic flow I f37l for 
different collision centralities within a hybrid hydro+cascade approach <38l and found that, at least with 
the standard Glauber model initial conditions, hadronic dissipation in the late hadron resonance gas stage 
can fully explain the deviations between data and ideal fluid dynamics at forward rapidities and in pe¬ 
ripheral collisions. Some additional QGP viscosity may be needed if the initial state is instead controlled 
by gluon saturation (see <38l for details). 

5.3 Where is the phase transition signature? 

The large hadron gas viscosity spoils one of the clearest experimental signatures for the quark-hadron 
phase transition, the predicted Q non-monotonic beam energy dependence of V 2 which was already 
described in Sec. Q and is shown in Fig. Eh. As one comes down from infinite beam energy, V 2 is 
predicted to first decrease (due to the softening of the EOS in the phase transition region) and then re¬ 
cover somewhat in the moderately stiff hadron gas phase. The hadron gas viscosity spoils this recovery, 


Pb+Pb, b=7fm 



Fig. 7: Left: Excitation function of radial and p ±-integrated elliptic flow for b = 7fm Pb+Pb or Au+Au collisions 
171 . The horizontal axis gives dN 7r /dy(b=7 fm), and the horizontal arrows indicate how dN^/dy was expected 
to correspond to the beam energy ranges covered by SPS, RHIC and LHC before the RHIC data were available. 
Right: Elliptic flow at fixed p± = 0.65 and 1.75 GeV/c from A+A collisions with A w 200 at a variety of beam 
energies <408 . 

leading to an apparently monotonous decrease of V 2 with falling beam energy (Fig.|fJ. However, recent 
PHENIX data from Au+Au collisions at y/s = 62 A GeV <40l indicate that this decrease may not be quite 
as monotonous as suggested by Fig. |6ji. Figure ITJr shows that, at fixed p±, the elliptic flow +2 (p ±) is 
essentially constant over the entire energy range explored at RHIC (from 62 to 200 A GeV), decreasing 
only when going further down to SPS and AGS energies. When integrated over p±, this turns into a 
monotonous behaviour as in Fig. due to the steepening p± spectra at lower y/s. While Fig. does 
not confirm the hydrodynamically predicted rise of V 2 towards lower y/s, it may still reflect this predicted 
non-monotonic structure in the elliptic flow excitation function, after strong dilution by hadronic viscos¬ 
ity effects [38 39:]. Obviously, many and more systematic hybrid calculations of the type pioneered by 
Teaney El are necessary to explore to what extent we can eventually prove the existence of a QCD 
phase transition using elliptic flow data in the SPS-RHIC energy domain. 











Obviously, it will be important to confirm non-viscous fluid behaviour at higher initial energy 
densities than so far explored, by extending Fig. |6jr to the right and verifying that v%/e settles on the 
hydrodynamic curve. This can be done with Pb+Pb collisions at the LF1C, or with full-overlap U+U 
collisions at RHIC 0. In addition, the large spatial source deformations achievable in full-overlap side- 
on-side U+U collisions allow for decisive systematic studies of the non-linear - path-length dependence 
of QCD radiative energy loss of fast partons. In 0 we give quantitative arguments why a U+U collision 
program should be seriously considered at RF1IC. 

5.4 Reconstructing the total momentum anisotropy from hadron data 

As discussed above in Section l5~2l no purely hydrodynamic model is able to simultaneously reproduce 
all low-px RHIC data, and (at the very least) some hadronic viscosity, implemented through a realistic 
hadronic cascade for the late expansion stage, must be considered. But if the observed hadron spectra 
and elliptic flow depend on the details of the late hadronic evolution ED, doesn’t this invalidate the claim 
that the momentum anisotropy is created early and probes the EOS of the hot QGP stage? 

The answer to this critical question is: “No, but_” Once the initial spatial eccentricity is gone, 

there is no further driving force for additional momentum anisotropy, and the only thing that can happen 
henceforth is that the latter is redistributed among the various hadron species and in transverse momen¬ 
tum (4dl . If the total momentum anisotropy saturates before the system hadronizes, hadronic viscosity 
is no problem for the overall momentum anisotropy (which is then entirely controlled by the preceding 
QGP and its EOS), and it can be fully reconstructed from the measured hadron spectra. However, the 
flow measure to look at for this purpose is not vi, but instead the yr -weighted elliptic flow (42). Using 
the kinetic theory relation T lll '(x) = f (d 3 p/E)p^ l p u fi(x,p) we can express the total final momentum 
anisotropy as a sum over contributions from all hadron species: 

fl „, = (T**-r»»)s, = [A M ‘ Ppi B> 

P {T~+T«) E| EMon J p 2 _A_ d 2 pi 

Here (.. ,)s f denotes an integral over the final kinetic decoupling hypersurface. Note that this expression 
is not entirely model-independent since it requires correcting the measured hadron spectra for post- 
freezeout resonance decays. 

For ideal fluids which live long enough for the spatial eccentricity to vanish and hence the mo¬ 
mentum anisotropy to saturate before the hadrons decouple, this e® nal is a clean probe of the (time- and 
temperature-averaged) speed of sound of the dense matter, as discussed in Sectional However, if the 
QGP hadronizes before e p saturates (i.e. if the system enters the hadronic stage with significant spatial 
eccentricity left over), ejj" al becomes sensitive to large dissipative effects during the late hadronic stage. 
This dilutes and distorts its sensitivity to the softening of the EOS (dropping speed of sound) during the 
quark-hadron phase transition, as discussed in the preceding subsection. An assessment of the remaining 
sensitivity of the excitation function (i.e. beam energy dependence) of 6p nal to the quark-hadron transi¬ 
tion and its usefulness as a phase transition signature |7)[39) thus requires systematic quantitative studies 
of the time evolution of e p within hydro+cascade hybrid models such as those used in 021 ESI- 

6 Ideal fluid response of the QGP to penetrating hard partons 

If the QGP indeed behaves as an almost perfect fluid, an interesting issue are possible hydrodynamic 
effects caused by jet quenching, i.e. by the localized energy deposited by a fast parton created early 
in the collision as it propagates through the dense fireball, colliding with the plasma constituents and 
radiating gluons E3- Due to space limitations I can only give a very abbreviated discussion here of 
this exciting question. There is experimental evidence (44) that the energy lost by fast partons traveling 
through the dense QCD medium formed in Au+Au collisions at RHIC largely thermalizes. Since the 
source of this deposited energy, the fast parton, travels at supersonic speed, it was suggested that this 
should generate either a hydrodynamic conical Mach shock wave propagating through the ideal QGP 
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fluid 1451 or a conical colored wake field propagating through the quark-gluon plasma OH S3- In either 
case, it was expected [45 46, |47|| that this should lead to anisotropic particle emission along a cone 
around the direction of the fast parton whose opening angle would reflect the speed of sound or the speed 
of collective plasma waves (both important plasma properties) in the quark-gluon plasma. 

A lot of excitement was generated by the fact that the PHENIX Collaboration l48l saw structures in 
the angular correlations of hadron emission relative to the direction of a fast trigger particle which might 
be evidence for such conical flow. The STAR Collaboration, on the other hand, did not see clear conical 
structures, but only a general broadening of the peak associated with the direction of the quenching jet 
144 491. A recent hydrodynamical simulation fl50l showed that conical flow generated by the fast parton 
and superimposed on the radial expansion flow of the fluid is indeed visible in the calculation, but that 
(even under the optimized conditions studied in E3) it does not manifest itself through the predicted 
peaks in the angular distribution of emitted hadrons. Other hydrodynamic effects, such as local heating 
and backsplash from the “crater” created by the fast parton, overlay the Mach cone phenomenon, the 
final result only being a large broadening of the peak associated with the quenching jet l50l . 

7 Conclusions 

The collective flow patterns observed at RHIC provide strong evidence for fast thermalization at less 
than 1 fm/c after impact and at energy densities more than an order of magnitude above the critical value 
for color deconfinement. The thus created thermalized QGP is a strongly coupled plasma which behaves 
like an almost ideal fluid. These features are first brought out in heavy-ion collisions at RHIC near 
midrapidity because only there the initial energy densities and QGP life times are large enough for the 
ideal fluid character of the QGP to really manifest itself, in the form of fully saturated hydrodynamic 
elliptic flow, undiluted by late non-equilibrium effects from the highly viscous hadron resonance gas 
which dominates the expansion at lower energies. 

We are now ready for a systematic experimental and theoretical program to quantitatively extract 
the EOS, thermalization time and transport properties of QGP and hot hadronic matter. This requires 
more statistics and a wider systematic range for soft hadron production data, but more importantly a 
wide range of systematic simulation studies with the “hydro-hadro” hybrid algorithms and, above all, a 
(3+l)-dimensional viscous relativistic hydrodynamic code (see E3 for more on that). 
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